Phase diagram of a probabilistic cellular automaton with 

three-site interactions 

A.P.F. Atman 0, Ronald Dickman I and J.G. Moreira i 
Departamento de Fisica, Instituto de Ciencias Exatas, Universidade Federal de Minas 

Gerais, C. P. 702 
30123-970, Belo Horizonte, MG - Brazil 
(February 1, 2008) 



Abstract 



We study a (1+1) dimensional probabilistic cellular automaton that 
is closely related to the Domany-Kinzel (DKCA), but in which the 
update of a given site depends on the state of three sites at the 
previous time step. Thus, compared with the DKCA, there is an 
additional parameter, p%, representing the probability for a site to 
be active at time t, given that its nearest neighbors and itself were 
active at time t — 1. We study phase transitions and critical behavior 
for the activity and for damage spreading, using one- and two-site 
mean-field approximations, and simulations, for P3 = and p% = 1. 
We find evidence for a line of tricritical points in the (pi, p2, Pz) 
parameter space, obtained using a mean-field approximation at pair 
level. To construct the phase diagram in simulations we employ 
the growth-exponent method in an interface representation. For 
Pz = 0, the phase diagram is similar to the DKCA, but the dam- 
age spreading transition exhibits a reentrant phase. For p 3 = 1, 
the growth-exponent method reproduces the two absorbing states, 
first and second-order phase transitions, bicritical point, and dam- 
age spreading transition recently identified by Bagnoli et al. [Phys. 
Rev. E63, 046116 (2001)]. 
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I. INTRODUCTION 



Probabilistic cellular automata (PCA) are widely used to model systems with 
local interactions in physics, chemistry, biology and social sciences Jl],[2|,||,[4]|| . 
Despite their simplicity, these models exhibit complex behavior and are used to 
investigate fundamental problems in statistical mechanics, such as spin models 
IK?] and nonequilibrium phenomena |||§. In particular, the problem of phase 
transitions in the presence of absorbing states has attracted increasing interest 



in recent years; PCA play a major role in these studies [jTQ|JTT] , [T^JT3|1 . The PCA 



introduced by Domany and Kinzel |§| is, along with the contact process [ fHHIS - 
one of the simplest models exhibiting an absorbing-state phase transition. 

The one- dimensional Domany-Kinzel stochastic cellular automaton (DKCA) 
is a completely discrete system - temporally, spatially and in its state space - 
which attracts interest as a particle system affording a test of ideas on scaling 
in nonequilibrium critical phenomena |16|]. The DKCA has a unique absorbing 
("vacuum") state; its phase diagram presents a critical line separating this ab- 
sorbing phase from an active one. Continuous phase transitions to an absorbing 
state are conjectured to belong generically to the directed percolation (DP) uni- 
versality class ||17|| . In addition to the active-absorbing transition, Martins et al. 
|18| found a damage spreading (DS) transition separating the active phase into 



nonchaotic and chaotic phases. There is numerical evidence that the critical be- 
havior along this transition line also belongs to the DP class, as expected on the 
basis of universality |19|| . 

Recently, Bagnoli et al || introduced a model that can be considered a natural 
extension of the DKCA: a one-dimensional PCA in which the update of a given 
site depends on the state of its nearest neighbors and itself, at the preceding 
time step. (We shall refer to this model as the BPCA.) Thus, compared with the 
DKCA, there is an additional parameter, p 3 , representing the probability for a 
site to be active at time t, given that all three sites were active at time t — 1. 
Bagnoli et al. studied p 3 = 1, in which case the model presents two absorbing 
states: the empty one and the completely occupied configuration. As in the 
DKCA, the density is the order parameter. These authors used the mean field 
approximation (at site level), simulations, and field-theoretic arguments to study 
the model. They found a rich phase diagram, with first- and second-order phase 
transitions, a bicritical point, and a damage-spreading transition. Except for the 
line P2 = 1 in the DKCA, this is the simplest PCA that exhibits a discontinuous 
phase transition M. 

In this work we extend the analysis of the BPCA considering two cases: the 
previously studied p 3 = 1, which corresponds to a ferromagnetic-like model, and 



p 3 = 0, representing a game- of-lif e-like model |20|j21||. We extend the mean-field 



analysis to the pair level, and use simulations to construct the phase diagram. In 
simulations, we apply the growth-exponent method [22] to identify transitions. 

This paper is structured as follows: in Section II we define the model and its 
interface representation; the site and pair mean-field approximations are discussed 
in Section III. Simulation results are presented in Section IV. We summarize our 
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findings in Section V. 



II. MODEL 

The one- dimensional PCA with three-site neighborhood (BPCA) was pro- 
posed by Bagnoli et al. |J. It consists of a ring of L sites (i = 1,2, ...,L), 
with periodic boundaries, in which each site % has two possible states, conve- 
niently denoted by <7j = 0,1. The state of the system at time t is given by 
the set {cjj(t)}. In contrast to the deterministic CA studied by Wolfram [|TJ, 
the present model is a discrete time Markov process: the rules for updating 
the system are given by transition probabilities. In particular, state of site i 
at time t + 1 depends on a^i(t), a^t) and o~m(t), via the transition probability 
P(<7i(t+l)|oi_i(t), (Ti(t), <7jfi(t)). The latter is of totalistic form, i.e., the depen- 
dence is through = Oi_i(t) +a#i(t). Since S(t) =0 implies cr f (t+1 ) = 
with probability 1, there remain three free parameters for defining the transition 
probability. Specifically: 
P(1|0,0,1) = P(1|0,1,0) = P(1|1,0,0) = Pl , 
P(1|0, 1, 1) = P(l|l, 0, 1) = P(l|l, 1, Q)=p2, 
P(l|l,l,l)=p 3 • 

Evidently, P(0\ai-i,a,a i+1 ) = 1 - P(l|(7 i _i, a i: a i+ x). 

Depending on the values of (pi,p 2 ,p 3 ), the asymptotic (t — > oo) state of the 
system is either in an absorbing phase (phase 0, with all sites in state 0, or 
phase 1, with all sites in state 1), or in the active phase, in which the stationary 
density p of sites in state 1 takes a value different from zero or one. Complete 
determination of the phase diagram in the three-dimensional parameter space is 
a rather difficult open problem. In this work we focus on two cases: p 3 — 1 and 
p 3 = 0. In the first case, the model possesses the two absorbing phases cited 
above, as well as an active phase and a chaotic region (associated with damage 
spreading). For p 3 = 0, the phase 1 is no longer absorbing (though the phase 
of course remains so), and there is again an active phase; the chaotic region is 
reentrant. p 3 = describes a situation in which "crowding" of individuals leads 
to their destruction, similar to Conway's Game of Life model PP],ZH|, while p 3 — 1 
corresponds to a ferromagnetic-like model. 

The absorbing/active transitions are continuous phase transitions, character- 
ized by critical exponents which belong to the DP universality class. The phase 
/phase 1 transition is discontinuous |J, and the exponents are those of compact 
directed percolation. The DS transitions are also in the DP class, consistent with 



Grassberger's prediction [19]. The termination of two critical lines at a line of 
discontinuous transitions marks a bicritical point, as has been found in the BPCA 
for p 3 = 1. For p 3 < 1, the phase 1 is no longer absorbing, so that one of the 
phase boundaries (i.e., between the active and phase 1-absorbing phases) is no 
longer present. We find (using the two-site mean-field approximation, discussed 
in Sec. IIIB) , that the bicritical point is actually one terminus of a line of tricrit- 
ical points: for each fixed p 3 in the range 1 > p 3 > p\ ~ 1/3, the absorbing-active 
transition is discontinuous for (p\,P2) < (p\,pi) and continuous for larger val- 



3 



ues, (pi,p2) > (Pi)P|)) where (p\,p\) is the tricritical point. For p 3 < p\, the 
absorbing/active transition is always continuous. 



A. Surface Representation 

The mapping of dynamical systems to a surface-growth representation is an 
interesting problem, since in many cases the resulting scaling properties are un- 
known a priori. Integration of the local activity (with respect to time) is the 
most natural procedure. The present method employs the interface representa- 



tion proposed by de Salles et al. p4J . The procedure consists in transforming the 
spatiotemporal patterns generated by the PCA to a solid-on-solid (SOS) particle 
deposition. The surface-growth process is attended by kinetic roughening; the 



associated critical exponents can be measured [II] following the scaling concepts 



developed by Family and Vicsek po[ . Atman and Moreira |22[ demonstrated 
that the growth exponent f3 w exhibits a cusp at criticality, and is very useful for 
detecting phase transitions. 

Height variables are defined by summing the variables &i(r) over the first t 
time steps: 

^(*) = I>(t) • (1) 

T=0 

In this way we generate a growth process, with correlations embodied in the 
roughness w(L,t) P5[| , defined by 



w 2 



(L,t) = ±(j:(h l (t)-h(t))^ , (2) 



where h(t) is the mean value of hi(t) at time t, and < . . . > denotes an average 
over realizations. 



For surface growth models [26], we expect w(L,t) to follow the scaling form 



w(L,t) ~ L a f , (3) 

where f(u) is a universal scaling function, a is the roughness exponent, z = 
a/ (5 W the dynamic exponent and j3 w the growth exponent. The function f(u) = 
constant for large u, while f(u) ~ for small u (t <C L z ). At short times, 
therefore, we expect w(t) ~ t^ . It is possible to measure (3 W from the slope of 
a log — log plot of w(L,t) versus t. In the active phase, the roughness does not 
saturate, growing instead as w(L, t) ~ t 1//2 , corresponding to uncorrelated growth 



22], since the correlation length is finite, away from the critical point. 



In previous work, Atman and Moreira [22] showed that f3 w attains a maximum 
at the phase transition, and measured its value along the transition line of the 
DKCA. DP values for others critical exponents in the surface representation of 
DKCA are verified in Ref. ??. 
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B. Damage Spreading 



Martins et al. used the damage-spreading technique to show that the 



active phase of the DKCA consists of two phases, chaotic and nonchaotic. The 
order parameter of this transition is the difference between two replicas initialized 
with different configurations, but subject to the same sequence of random events 
during the subsequent evolution. In practice, we let the system evolve until it 
attains a stationary state, and then copy the configuration, introducing some 
alterations ("damage"). The two replicas, one with state <Ji{t) and the other 
with state Qi(t), then evolve using the same sequence of random numbers, and 
the difference between them, 

T l (t) = \<Ti(t) - Qi(t)\ , 

is monitored. The fraction of sites in the two replicas with <7j 7^ Qi defines their 
Hamming distance: 

D H (t) = . (4) 

i 

The stationary Hamming distance is null in the non-chaotic phase and positive 
in the chaotic phase. 

To study the chaotic/non-chaotic boundary, we use a slightly different method, 
in which the difference between the replicas is used to generate a surface growth 
process, as described above. In this case, 

k(f) = Er,(T). (5) 

T=0 

Thus, the profile generated by the difference between the replicas behaves just as 
the profiles at the phase 0/active boundary: the roughness reaches a stationary 
value in the non-chaotic phase and grows indefinitely in the chaotic phase, with 
a cusp in the j3 w value at the transition. 

Since the system has already relaxed to the stationary state when we initiate 
the damage experiment, creating the damage by randomly altering sites in one 
copy is likely to perturb the particle density and correlation functions away from 
their stationary values. This in turn would introduce an undesirable assymetry 
between the replicas, since the dynamics of damage spreading would be mixed 
with that of relaxation back to the stationary state (in the copy but not in the 
original). To avoid such complications, we generate the damage by rotating the 
copy by 180° with respect to the original, with no further modifications, that is, 
g{i,to) = a(i + L/2,to), subject to periodic boundary conditions. This represents 
a large initial damage (a Hamming distance of ~ 2p(l— p), with p the stationary 
particle density), which is statistically uniform over the system. 
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III. MEAN-FIELD THEORY 



A. One-site approximation 

To begin, we present the mean-field approximation at the site level, for p 3 = 
1 and P3 = 0, and construct the phase diagram from the equations obtained 
in this approximation. The BPCA is a Markov process in which all sites are 
updated simultaneously. The configuration {cr} is a set of stochastic variables 
with probability distribution at time t given by Pt{o). The evolution of the latter 
is governed by 

Pv rl =Y,u(*W)Pt{p>) , (6) 

a' 

where u(a\a') denotes the probability of the transition a' — > a, with the properties 
uj(a\a') > and J2a u(a\a') = 1. The transition probability for the BPCA is a 
product of factors associated with each site: 

L 

uj(a\a') = Y[wi(ai\a') , (7) 

i=l 

where Wi(ai\a') > is the conditional probability for site % to be active at time 
t + 1 given the configuration a' at the preceding step. The probabilities Wi are 
translationally invariant and in fact depend only on the states <Jj_i, Oi and a i+ \ 
at the previous step: 

Wi(ai\a') = wssfaWi^, v'i, . (8) 

We list the w 3s in Table I. 

Of interest are the n-site marginal probabilities. The evolution of the one-site 
distribution Pt(&i), is given by 

Pt+i{cri) = Him W3 s (o-ik-i, cr-, rt+JPtW-!, cr'i, cr' i+1 ) , (9) 

a'. , a 1 , a'. 

2—1 % l+l 

where Pt^cr'^, a'^ a' i+l ) is the marginal distribution for a set of three nearest- 
neighbor sites. The evolution of the two-site distribution is given by: 

Pt+i(cri,cr i+1 ) = J2J2J2J2 ^(^K-i^^'+i) 

a 'i-l a 'i a 'i+l a 'i+2 

X w 3a (a i+1 \<T , i , cr' i+1 , <^ +2 )P t {o , i _ 1 , cr-, cr' i+1 , a' i+2 ) . (10) 

Evidently we have an infinite hierarchy of equations. In the n-site approximation 
the hierarchy is truncated by estimating the (n+l)-site (and higher) probabilities 
on the basis of those for n sites. 

The simplest case is the one-site approximation, in which P t (cr^_ 1 , o[, cr' i+1 ) is 
factored so: P t (<7^_ 1 , cr-, <r' i+1 ) — P*(cr|_ 1 ) Pt(cr-) Pt(cr- +1 ) . This yields the recurrence 
relation, 
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Pt+1 = P3pt + 3P2pt(l - Pt) + 3pi/3t(l - ptf , 



(11) 



where p t = -Pt(l) is the density of active sites (i.e., the order parameter). 

Depending on the value of (pi, p%, p 3 ), Eq. ([□]) admits different stationary- 
solutions, corresponding to the possible BPCA phases discussed above: phase 
(p = 0), phase 1 (p = 1) and active (0 < p < 1). In order to verify the stability of 
the stationary solutions, we consider a small perturbation in the stationary value 



p* , p t = p* + Ap t . Applying this variable change in the Eq. (JTT) , we obtain for 
the mean-field approximation at the site level, 

Ap t+1 = A Pt [3p* 2 (p 3 + 3 Pl - 3p 2 ) + 6p*(p 2 - 2 Pl ) + 3 Pl ] + 
(A Pt ) 2 [3p*(p 3 + 3pi - 3p 2 ) + 3(p 2 - 2 Pl )] + 
(Api) 3 [p3 + 3p!-3p 2 ] . (12) 

We can write Eq. ([TJ) in a simplified manner, 

Ap t+1 = a(p*)Ap t + b(p*)A 2 p t + cA 3 p t , 

where the coefficients a(p*), b(p*) and c can be associated with the stability of 
the solutions. 

Considering the solution, p* = 0, the stability condition is a(p*) < 1, which 
implies p\ < 1/3. In the case a(p*) = 1, (pi = 1/3), the solution will be stable 
only if 6(p*)<0 and c<0; the first condition implies p3<2/3. For a(p*) = l and 
b(p*) = 0, the stability the condition c < leads to p 3 < 1. Thus, the solution 
p* = is always stable for pi < 1/3, p 2 < 2/3 and ps < 1. The point (1/3, 2/3, 1) 
corresponds to a tricritical point in this aproximation. 

Considering the solution, p* = 1, the stability condition a(p*) < 1 implies 
P3 — T>2 < 1/3. In case a(p*) = 1, (p 3 =p 2 + 1/3), the condition b(p*) < yields 
Vi ~ Pi > 1/3. For a(p*) = 1 and b(p*) = 0, the stability condition c < implies 
P3<1. So, the point (1/3, 2/3, 1) is also a tricritical point for the solution p* = l! 
Thus, the point (1/3, 2/3, 1) corresponds to a bicritical point, as Bagnoli et al. 
have already shown. In fact, the solution p* = 1 is absorbing only for p 3 = 1 (since 
for p 3 < 1, the dynamics of updating destroys the phase 1), and this solution is 
stable for p 2 > 2/3. 

For < p* < 1, the stability condition a(p*) < 1 implies a inequality of 
second degree in terms of p*. Considering a = 1, we can solve the corresponding 
equation, and obtain 



(2pi -p 2 )±J (p 2 - 2pi) 2 - (p 3 + 3pi - 3p 2 )(3pi - 1) 

p = . 

P3 - 3pi - 3p 2 

Note that p* vanishes on the line (1/3, 2/3, P3). It is easy to see that the positive 
solution, p* +) is valid for any (pi, p 2 , p 3 ), but the negative solution, pi, is valid 
only for p 2 < 2p\. Considering the plane p\ = 1/3, the transition line for p 2 > 2/3 
coincides with the vanishing of the square root, since p* must be real. This implies 
a discontinuous transition for p 2 > 2/3, and the line (1/3, 2/3, p 3 ) corresponds to 
the tricritical line in (p 1; p 2 , p 3 ) space, in the site approximation. 
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In the simulations we consider two cases: p% = 1 and p 3 = 0. Considering the 
stability analysis above, we can summarize the phase diagram in these two cases 
in the following way: for p 3 = 1, Eq. ( O ) can be written as, 

P ((3pi - 3p 2 + l)p 2 + (3p 2 - 6 Pl )p + 3 Pl - l) = . (13) 
The three solutions of this equation are: 

• p = - Phase 0, stable for p% < 1/3 and p 2 < 2/3; 

• p = 1 - Phase 1, stable for p 2 > 2/3 and p\ > 1/3; 

• Active phase, for p 2 < 2/3 and p\ > 1/3, where the stationary density is 
given by: 

P 3pi - 3p 2 + 1 ^ 

For p 2 > 2/3 and p\ < 1/3, we have a discontinuous transition line separating 
the phase 0/phase 1 (both stable) at p 2 = 1 — p%. 
In case p 3 = 0, we have: 

P ((3pi - 3p 2 )p 2 + (3p 2 - 6p x )p + (3 Pl - 1)) = . (15) 

For this equation there are only two distinct phases: 

• p = - Frozen phase, stable for pi < 1/3, p 2 < 2/3; 

• Active phase, valid for p\ > 1/3 and pi 7^ p 2 , where the stationary density 
is given by: 



6pi - 3p 2 ± ppi - 12p 2 + 12 Pl 

P = R R ^ 16 

6pt - 6p 2 

the negative root is valid for p 2 < 2p±, while the positive root is valid for 
p 2 > 2pi. 

For p 2 > 2/3, the solution of eq. ([16]) is either complex or strictly > 0, implying 
a discontinuous transition from the phase to the active phase, as we antecipate. 

The phase diagram in the one-site approximation, for p 3 = 1, is shown in 
Fig. 1; that for p 3 = is shown in Fig. 2. For p 3 = 1, the phase diagram is as 
expected [[J. For p 3 = 0, we expect a single absorbing state, but the behavior for 
p 2 > 2/3, where the transition line is discontinuous is not expected on the basis 
of simulations. (At the transition the r.h.s. of Eq. (O) changes from a complex 
to real, nonzero value.) 
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B. Pair approximation 



At the pair level, the probability P t ((T^_ 1 , a[, cr' i+1 , * s factored in the fol- 
lowing way 

Calling z = P(l,l), and noting that P(l) = P(1,0) + P(l, 1), we can write 
P(l, 0) = P(0, 1) = p — z, and P(0, 0) = 1 — 2p + z. The recursion relations for 
the density of active sites p and for the density of active pairs z are: 

Pt \Pt 1 - Pt J 

( 2(1 -2p t + zt) ( P t-zt) \ 
+ Pl(pt-Zt)[ : + , 18 

V i - pt pt J 



Zt+i = v\ 



'(p t -z t ) 2 (l-2p t + z t )(2-p t y 



ptO- - pt) 



2 



+ 2 PlP2 ( ^ Zt \ ) (z t (l - 2p t + z t ) + (p t - z t ) 2 



,Pt( l ~ Pt), 

+ ft I rf(i- ft ) j + 2KK l^^J + ' (19) 

Iterating these relations numerically until a steady state is reached, we construct 
the phase diagram in the pair approximation. Results for p 3 — 1 and p% = are 
shown in Figs. 1 and 2, respectively. 

Using the pair approximation, Eq. (|18D and Eq. (0), we find numerically the 
critical surface in the (pi,P2,P3) parameter space and the line of tricritical points, 
for p 3 < 1, as sketched in Fig. 3. In the region p 3 < 1, the phase 1 disappears and 
the discontinuous absorbing transition lines meet the continuous transition lines 
at the tricritical points, as shown in Fig. 3. Each one of the eight vertices in this 
diagram corresponds to a different deterministic rule in the automata studied by 
Wolfram [[!]]; for example, (p\ = 1, p 2 = 0, p 3 = 0) corresponds to the rule 22; 
(pi — 1, P2 — 1, P3 = 0) to rule 126; etc. 



C. Damage Spreading Transition at Site Level 

Bagnoli et al. || derived a mean-field approximation for the DS transition at 
Ps = 1, showing that there is a chaotic region in the active phase of the BPCA. To 
obtain the one-site approximation for the BPCA at p 3 = 0, we use the approach 
of Tome [p3|j ; denoting the configurations by {a{\ and {t^}, the Hamming distance 
is given by 

H t = ({a l -r l f) , (20) 
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where the brakets denote an average over realizations. The evolution of the joint 
probability follows 

P m (a; r) = £ W(a; r \ a'; r')P t (a'; r') , (21) 

where 

W(a; r | a'- r') = J[ wfa r, t \ o[_ x , a[, a' i+1 ; r^, r(, r' l+1 ) , (22) 

i 

is the transition probability for two the two systems (subject to the same noise), 
from the state (a'; t') to (a; r). Using the transition probabilities defined in Table 
I, we can calculate the joint transition probabilities, as shown in Table II. 

Now, we can write the equations for the evolution of the order parameter 
associated with the chaotic transition - the Hamming distance. Denoting the 
Hamming distance defined by Eq. (|J) as if> t = P*(l; 0) = Pt(0; 1), and using the 
relation (^T|), we have 

ipt+i =< E7(l;0 | (7j_i,(7i,cr i+ i;Tj + i,Tj,Tj + i) > ; 

using the rules of Table II, we can write if> t as 

P m (l; 0) = p 3 P(l, 1, 1; 0, 0, 0) + 3p 2 P t (l, 1, 0; 0, 0, 0) + (23) 
+ 3(6 + c)P t (l, 1, 1; 1, 1, 0) + 3 Pl P t (l, 0, 0; 0, 0, 0) + 
+ 3(6' + c')P(l, 1, 1; 1, 0, 0) + 9(6" + c")P t {\, 1, 0; 1, 0, 0) . 

Calling -Pt(l) = x t , we can write P*(l; 1) = x t — if>t and Pt(0; 0) = 1 — x t — if>t] 
thus, using the one-site mean-field approximation P t (<7j_i, <jj, 0i+i; Tj_i, Ti, r i+ i) = 
P t ((T i _i,(T i ,(T i+ i)Pt(T i _ 1 ,T i ,T i+ i), we can write Eq. (|3]) as 

1>t+i = MPstf + 3P2^(1 ~xt- A) + 3pi(l - x t - iptf + 3(6 + c){x t + 4j t ) 2 (24) 
+3(6' + c')ip t (x t - if) t ) + 9(6" + c")(x t - if; t )(l - x t - Vt)] • 

Finally, considering the case p 3 = 0, we have (6 + c) = p 2 , (6' + c') = pi and 
(6" + c") =| pi — P2 |; inserting these values in Eq. (|24| ) we obtain 

V>m = 3^{3 | px - p 2 | ^ 2 + [(p 2 - 2px) - 3 | p a - p 2 | +3(p x - p 2 )x^\) t (25) 
+ [(Pi +P2 - 3 | pi -p 2 \)x 2 t + (3 I pi -p 2 | -2px)x t +pi]} . 

This equation can be iterated numerically using the stationary values of Xt ob- 
tained from Eq. (0); there are three possibilities for the joint solutions of Eq. 
(jnj|) and Eq. (|25|): x = ?/> = 0, corresponding to the phase 0; x = 0, if) ^ 0, 
corresponding to the active phase; and x, ifi j£ 0, corresponding to the chaotic 
phase. In Fig. 4 we show the stationary solutions for these equations for some 
values of p 2 . We note that there is a discontinuous DS transition line in this 
approximation: for p 2 >pt (P2 = 0.9 in Fig. 4, for example), if> is always positive if 
x > 0. It implies that the DS transition line for p 2 >Pt falls on the discontinuous 
transition analysed in subsection III A. 
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IV. SIMULATION RESULTS 



We construct the BPCA phase diagram, for p 3 = and p 3 — 1, using simu- 
lations of systems of up to L = 10000 sites (with periodic boundaries), applying 
the growth exponent method \22\ to locate the transition lines. The initial con- 
dition used in the simulations is random, with half the sites occupied. The phase 
diagrams for absorbing-state transitions are shown in Figs. 1 and 2, for p 3 = 1 
and p 3 = 0, respectively. As expected, the pair approximation yields a better 
prediction than does one-site mean-field theory. Note that for p 3 = 0, the phase 
diagram is qualitatively the same as for the DKCA the major difference 

is that the active phase is quite enlarged in the BPCA. 

For ps — 1, the pair approximation prediction for the phase boundaries is 
qualitatively correct, although the bicritical point remains in the same position 
(pi = l/3,p2 = 2/3), as in the site approximation. Simulations place the bicritical 
point at (0.460(3), 0.540(3)), but the phase boundaries are in reasonable agree- 
ment with the pair approximation prediction. It is important to note that there 
are only three transitions in this diagram: the phase /phase 1 transition (discon- 
tinuous) and the phase 0/active and active/phase 1 transitions (continuous). All 
transitions are located using the growth-exponent method, confirming that this 
method is able to detect both continuous and discontinuous phase transitions. 

The phase diagrams for the DS transition are shown in Figs. 5 and 6, for p 3 = 
and p 3 — 1 respectively. In the case p 3 = 1, we confirm the results of Bagnoli 
et al, but some comments are in order. In Ref. f9j, the authors sketched several 
"damaged domains" that appear along the active/absorbing phase boundary, 
and attributed them to the divergence of the relaxation time, or to the fact that 
small differences in the initial configuration can drive the system to a different 
absorbing state. As shown in Fig. 7, where we compare the DS transitions 
obtained using (1) "rotation" damage and (2) random damage in 10% of the sites, 
these domains are only associated with the absorbing-state transition. We see two 
maxima in the (3 W x pi curves in Fig. 7: the left maximum, more pronounced, 
that corresponds to the absorbing/active transition, and the right maximum, 
which actually corresponds to the nonchaotic/chaotic transition and belongs to 
the DP class. The left maximum can yields to an "apparent" DS transition, as 
the 'islands' of damage commented in the Ref. ||. 

Thus, the left maximum corresponds to the "damage domains" of the Ref. 
||, but in fact it is due the absorbing/active transition: in this region, when the 
replica is created (in the stationary state), it turns out that only a small region 
of the ring is active, which results in a constant contribution to the Hamming 
distance (proportional to 2p* with rotation damage, and to 2p(t), with random 
damage). Thus, with rotation damage, the exponent (3 W approaches 1 because 
the system is already in the stationary state, while with the random damage, 
(3 W ~ 0.84, due decay of the activity, recovering the DP value. This behavior 
should not be confused with the true damage transition that occurs only at the 
second maximum, and which corresponds to a unique phase boundary, as shown 
in Fig. 6. 
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In the case p% = 0, the simulations confirm the reentrant chaotic transition 
predicted by the one-site mean-field approximation. As shown in Fig. 5, the DS 
transition line is concave inward toward the active phase, and presents distinct 
behaviors for p% > p 2 and p\ < p 2 , as expected. We note that the simulation 
results suggest a sudden cnahge in the orientation of the active/ chaotic phase 
boundary where the latter crosses the line p 2 = Pi- Such discontinuity of slope, 
clearly evident in the MF prediction, may be a consequence of the singular be- 
havior of several transition probabilities on p\ — p 2 parameters space, as shown 
in Table II and Eq. 

V. CONCLUSIONS 

In this work we apply the growth-exponent method in Monte Carlo simula- 
tions, and one- and two-site mean-field approximations, to construct the phase 
diagram of the BPCA for p% = 1 and p^ = 0. The method detects both first- and 
second-order phase transitions, and also can be used to locate DS transitions. 
The exponent values indicate that all continuous phase transitions belong to the 
directed percolation universality class, while the exponent at the discontinuous 
phase transition agrees with the compact directed percolation value. 

We find evidence of a line of tricritical points in the (pi, p 2 , P3) parameter 
space, using the mean-field pair approximations. We also find a reentrant chaotic 
transition for ^3 = in the mean field approximation, that was confirmed by 
simulations. These observations illustrate the rich and at times surprising phase 
space structure found in simple nonequilibrium systems. 

Acknowledgements We thank the Brazilian agencies CNPq and Fapemig for 
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TABLES 



TABLE I. BPCA transition probabilities 

cri/K^X+i) (1,1,1) (LLP), (LOT) or (0,1,1) (0,0,1), (L0,0) or (0,L0) (0,0,0) 

1 P3 P2 Pi 

1 - p 3 1 - p 2 1 - Pi 1 



TABLE II. Joint transition probabilities for two BPCA subjected to the same noise. 



Pi, n 


(LLLLLi) 


(LLLLL0) (i,i,i;i,0,0) 


Li 


P3 


min(p 2 , P3) min(pi,p 3 ) 


1,0 





b = max (p 2 - p 3 , 0) b' = max (p 3 - pi, 0) 


0,1 





c = max (p 3 - pi, 0) c' = max (p 1 - p 3 , 0) 


0,0 


1 ~P3 


l-max(p 2 ,P3) l-max(pi,p 3 ) 




Pi, n 


i,i,i;0,o,o 


(1,1,0;1,L0) (LL0;1,0,0) 


Li 





P2 min(pi,p 2 ) 


1,0 


P3 


6" = max (pi — p2, 0) 


0,1 





c" = max (p2 — pi , 0) 


0,0 


1 ~P3 


1 -p 2 1 - max(pi, p 2 ) 




Pi, n 


(1,1,0;0,0,0) 


(L0,0;1,0,0) (L0,0;0,0,0) (0,0,0;0,0,0) 


Li 





Pi 


1,0 


P2 


pi 


o,i 








0,0 


1-P2 


1 - Pi 1-Pl 1 
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Figure Captions 
Fig. 1. Phase diagram for the BPCA, p 3 = 1. One- and two-site mean- field 
approximations are compared with simulation results. 

Fig. 2. Phase diagram for the BPCA, p 3 = 0. The mean field approxima- 
tion predicts a tricritical point at (p 2 = 2/3, p\ = 1/3), where a discontinuous 
boundary (dashed line) meets a continuous transition line. 

Fig. 3. Line of tricritical points and the critical surface in the (pi, p 2 , P3) 
parameter space, as predicted by pair level mean field approximation. The tri- 
critical line ends near p 3 ~ 1/3, while at the site level it extends to p 3 = 0, as 
shown in Fig. 2. 

Fig. 4. Density of active sites and Hamming distance, in the mean-field ap- 
proximation, at site level, for p s = 0. Note that for p 2 > 2/3, the transition is 
discontinuous. 

Fig. 5. DS transition line, for p 3 = 0. The one-site mean-field approximation 
(inset) is compared with simulation data (main graph). Reentrant behavior is 
observed in both cases. 

Fig. 6. DS transition line for p% = 1. The DS boundary meets the absorbing 
transition lines for (p 1? p 2 ) near 1. 

Fig. 7. DS growth exponent for two different initial damages, for p 3 — l. The 
left maximum corresponds to the absorbing transition and the right maximum 
corresponds to the DS transition. 
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